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Abstract 

Residual wavefront errors in optical elements limit the performance of coronagraphs. 
To improve their efficiency, different types of devices have been proposed to correct 
or calibrate these errors. In this paper, we study one of these techniques proposed 
by Baudoz et al. 2006 [1] and called Self-Coherent Camera (SCC). The principle of 
this instrument is based on the lack of coherence between the stellar light and the 
planet that is searched for. After recalling the principle of the SCC, we simulate 
its performance under realistic conditions and compare it with the performance of 
differential imaging. 
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Etude des Performances d'une Self-Coherent Camera. La qualite de surface 
des optiques limite les performances des coronographes. La correction ou I'etalonnage 
de ces defauts optiques permet d'en ameliorer I'efficacite. Nous etudions dans cet 
article une technique, proposee par Baudoz et al. 2006 [1], qui permet d'etalonner 
les tavelures creees au plan focal par les defauts de surface d'onde. Le principe de cet 
instrument, appele Self-Coherent Camera (SCC), est base sur I'absence de coherence 
entre I'etoile et la planete. Apres un rappel du principe de la SCC, nous presentons 
une comparaison de ses performances avec celles de I'imagerie differentielle. 
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1 Introduction 



Radial velocimetry has enabled the indirect detection of more than 200 ex- 
oplanets over the past decade [2]. The study of their physical parameters 
involves the measure of their spectra and a straightforward solution is their 
direct detection. However, even the brightest planets [3] are 10"^ to 10^ times 
fainter than their host and often located within a fraction of an arcsecond of 
the star. A large number of coronagraphs have been proposed to suppress the 
host star's overwhelming flux [4] but all of them are limited by the imper- 
fections of the wavefront. Using high-order Adaptive Optics (AO) or spatial 
observatory, the performance of coronagraphs is still limited by the aberrations 
in the coronagraphic optics and the residual errors from the AO (for ground- 
based observation). Thus, when attempting to detect a faint companion on 
long exposure images, residual speckle patterns can be a dominant source 
of error [5]. Since these speckle patterns are slowly drifting, it is mandatory 
to find a way to discriminate the speckles of the star from a faint compan- 
ion during the exposure. This is the purpose of differential imaging techniques. 
Several criteria have been proposed so far to discriminate the speckles from the 
planets: spectrophotometry [5,6], polarimetry [7,8], and coherence [1,9,10,11]. 
While both concepts based on spectrophotometry and polarimetry depends on 
the physical properties of the planets, the coherence is a robust criterion when 
no physical information is available from the planets that could be observed. 
Here, we study the way to calibrate the speckles using coherence as proposed 
by Baudoz et al. 2006 [1]. First, we recall the concept of the technique, called 
Self Coherent Camera (SCC). Then, we compare the detectability of a com- 
panion using the SCC with classical techniques of differential imaging. 



2 Principle 

The concept of the SCC has already been presented [1] but the main features 
are recalled here. The purpose of that device is to discriminate in a field 
of view a companion image from the speckles created by wavefront defects. 
While both features look almost the same on the detector, only the speckles 
are coherent with the stellar beam. A Fizeau recombination is used to encode 
the field of view with a coherent fringe pattern that affects only the stellar 
speckles. The principle and a possible set-up are described in figure 1. The 
light coming from the telescope is split into two beams. One of the beams is 
spatially filtered using a pinhole or an optical fiber. The typical size of the 
pinhole is about the size of the Point Spread Function core (^ ^/D) for two 
reasons: 1) Since high frequency pupil defects are diluted in the focal plane, 
the pinhole cleans the wavefront of the reference beam [12]. 2) The companion 
is not transmitted whenever its distance to the star is larger than the size 
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of the pinhole. The two beams are recombined in the focal plane in a Fizeau 
scheme. To do so, the two pupil beams are optically brought at the same plane 
right before a focusing lens. Because of the Fizeau interferences, the intensity 
distribution of the stellar flux is fringed. Because of the spatial filtering, the 
pupil illumination of the reference beam is not anymore uniform. Thus, the 
fringe contrast will not be 100%. The contrast will also be diminished by 
differential aberrations between the reference and the main beam. However, 
section 3 and 4 show that it does not limit the SCC performance. Since the 
fiux of the companion is removed from the reference beam by spatial filtering, 
the intensity of the companion will be unaltered by the reference beam and 
the image of the companion will not be fringed (figure 1). 



3 Formalism 

The electromagnetic field in the entrance pupil plane is described by ^1^(0 = 
A{^). The term A{x) is the complex amplitude in the focal plane and A{^) 
denotes the Fourier Transform of A{x). Coordinates ^ and x are used for the 
pupil plane and the focal plane respectively. The field of both the star and its 
companion in the pupil plane can be written :^'*(^) + ^<7(^) = -A*(^) +Ac{^). 
The coherent pupil of the reference beam can be described in the pupil plane 
by: ^niO = MO 

Assuming that the vector describes the distance between the coherent pupil 
and the corrugated pupil (supposed to be centred on zero), the field at the 
SCC pupil output plane is given by: 

*5cc(0 = **(0 + *c(0 + *i?(0 * s{c - Co) (1) 

where * is the convolution symbol. For a given complex amplitude A{x), we 
write the intensity detected at the focal plane as I{x) —\ A{x) p. Assuming 
the optical path difference of the interferometer is kept at zero, the image 
recorded at the focal plane is: 

I{x) =1 r= + Ic{x) + Ir{x) + 2Real{A4x)A*j^{x)} ■ cos(27ra;^o) (2) 

The image looks like a classical Fizeau image with fringes pinning the intensity 
distribution of the main star. However, the image of the companion does not 
show any fringes because it is not coherent with the reference beam. Since 
the intensity of the main star is coded with fringes while the image of the 
companion is not, it looks clearly that one can discriminate the image of the 
companion from the stellar flux. The simplest approach is to use the FFT 
of the image to separate the three areas limited by the autocorrelation of the 
pupil function. The centred area is the sum of the autocorrelations of the three 
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Figure 1. Left: Possible Set-up for the SCC. Right: Principle of the SCC. The elec- 
tromagnetic field distribution in the SCC pupil is indicated for the star and its 
reference beam (black) and a companion (grey). At the bottom right, the numerical 
Fourier Transform of the detected image shows autocorrelations and intercorrela- 
tions between the reference beam and the stellar beam. 

images Ir, h and Ic- The correlation between the fields and ^ r appears 
in lateral areas (Fig. 1). The separation between the two images of the pupil 
^0 is large enough so that the different terms are not superimposed and can 
be numerically separated. 

Assuming that we describe the centred area by Ii = 1^ + Ic + Ir and the two 
other areas by /_ = A*^.A^ and J+ = (/-)*, we can write the intensity of the 

companion using the following formula Iq = h — Ir . 

Ir 

To detect the companion Ic, one needs to record separately Ir. It can be done 
because Ir is a spatially filtered beam that can be very stable over time. In fact, 
the stability of Ir is a possible limitation of the SCC. A laboratory experiment 
is under development to analyze the impact of Ir instability. Assuming we 
record Ir, which is an estimate of Ir, the final Ic estimator, called SCC, is 
given by : 

SCC = h-hi-h^ (3) 

This equation is valid only for an exposure shorter than the coherent time of 
the atmospheric speckles and for monochromatic beam [1]. However, using a 
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chromatic corrector as proposed by Wynne 1979 [13]. the bandwidth could be 
increased up to a reasonable value (R=5). Equation 3 is true for any aberration 
or illumination of or 1^. Thus, neither the non-uniform illumination of Ir 
because of spatial filtering nor any differential aberration between the reference 
and the main beam should limit the detection of a companion with the SCC. 



4 Numerical simulations 



To confirm that result, we developed a numerical code that simulates the SCC. 
In this section, we decribe the hypothesis we assumed and detail the data 
processing of the simulation. The main physical hypothesis and parameters 
are listed underneath : 



• The telescope is 8-meter diameter. 

• The considered star is a 5th visible magnitude G star observed at 0.8//m 
with = ^ = 8 assuming a perfect Wynne corrector is used[13]. 

• The quantum efficiency of the camera is set to 0.4. The exposure time must 
be shorter than the coherence time, so we choosed an exposure time of 6 ms. 
As shown by Sarazin and Tokovinin's study [14], the atmospheric coherence 
time at the wavelength of 0.8 iim is longer than 6 ms more than 50% of the 
time at Paranal. 

• We used the approach proposed by Rigaut et al. 1998 [15,16] to generate the 
atmospheric phase screens of the Paranal site. Shack-Hartmann wavefront 
sensor with 80 actuators across the telescope diameter and a 1 ms AO closed 
loop temporal have been assumed. The seeing has been set to 0.6 arcsec at 
0.8 ^m. A global AO Power Spectral Density (PSD) taking into account 
servo-lag, aliasing and fitting errors has been computed. From this DSP, 
independent phase screens realizations are created. 

• Following Cavarroc et al. 2006 [17], we introduced common 5c and non- 
common 5nc static aberrations. For the simulated SCC (respec. differential 
imaging) device, 5c are static aberrations in the instrument upstream of 
the SCC (respec. differential imaging instrument). 5]sic are the differential 
static phase aberrations between the corrugated and the reference beams 
for both techniques. 5c is set to lOnm rms and 5]sic to 5nm or Inm rms 
and both are created from a PSD following a law [18]. 

• Photon noise, but no read-out-noise, is taken into consideration. 

• Images are 512x512 pixels and the entrance pupil diameter is D = 80 pixels. 
The simulation is monochromatic. 
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4-1 Self- Coherent Camera 



To simulate the SCC we split the telescope AO beam into two beams (see 
section 2) which are both corrugated by the AO corrected aberrations and Sc- 
The reference beam is filtered by a diameter pinhole in a focal plane. 
Then, an entrance pupil size diaphragm is used to stop the diffracted light 
into the following pupil plane. Finally, that reference pupil is corrugated by a 
non-common static phase aberrations, S^c- Then, the two beams, separated 
by 2.15 D, are recombined in a Fizeau scheme assuming a zero optical path 
difference. The resulting focal plane intensity is recorded in a numerical image. 

The SCC data processing is done in three steps. The first one is a Fast Fourier 
Transform of the interferometric image. Then, the three autocorrelation ar- 
eas, /i, /+ and /_, are separated using masks and an inverse Fast Fourier 
Transform gives /i, /+ and /_. Finally, the SCC residual image is given by 
equations, where J/j is a long exposure recording of the reference image. How- 
ever, the division by Ir leads to undefined values where Ir equals zero. Thus, 
the image SCC.Ir is computed instead of the SCC one. Afterwards, we divide 
the resulting image by Ir using a 10~^ threshold. 



4.2 Differential imaging device 

Both beams of that device are corrugated by the AO corrected aberrations 
and 5c- One of the beams, called here the reference beam for SCC comparison, 
is also corrugated by S^c- The differential image is the subtraction of the 
reference image to the other one. The simulation of differential imaging is 
optimistic because both beams are supposed observed at the same wavelength 
and at the same polarization. Moreover, the companion is introduced in only 
one image and completely removed from the other beam. 



4.3 Detectability 

The goal of the two previous simulated devices is to detect a faint compan- 
ion. Calling C the companion energy contrast, the image maximum intensity 
is C hmax, where hmax is the host star image maximum intensity. To deter- 
mine which contrast can be detected by the considered device, we define D, 
the detectability at 5a, as D = 5 . ^ , where N is the residual noise into the 
final image when no companion is present. In that paper, is computed as 
the square root of the azimutal spatial variance of the final image. 
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Figure 2. SCC (solid line) and differential imaging (dot-dashed line) detectabilities 
versus angular separation . Theoretical photon noise (dotted line) and 5nm rms 
non-common static aberrations (dashed line) are also plotted. The SCC reference 
image exposure time is 18 s. Left: One single image ofQms. Right: 10000 images 
of 6ms. 



5 Results of the simulation 



5. 1 Detectability versus angular separation 



Figure 2 shows SCC (solid line) and differential imaging (dot-dashed line) 
detectabilities at 5a versus the angular separation to the host star. For the 
SCC, the J/j exposure time is 18 s to minimize the photon noise. Using a single 
image (Figure 2 left), the theoretical photon noise (dotted line) is greater than 
the 5nm rms non-common static aberrations (dashed line). The SCC and 
differential imaging detectabihties roughly follow the theoretical photon noise 
limit. However, the SCC profile shows spikes and an amplification for angular 
separation greater than 8-^. These spikes and amplification are due to the Iji 
division where Iji is smaller than one photon. For a sum of 10000 images, the 
theoretical photon noise is decreased by a factor 100 (figure 2, right). Thus, 
the non-common static aberrations start to dominate in the differential image 
as predicted by Cavarroc et al. 2006 [17] . The SCC noise almost follows the 
photon noise limitation and is not limited by static aberrations where the 
reference image J/j is greater than one photon. However, spikes are greater 
than on the single image profile (figure 2, left). Other data processing are 
under study to minimize these spikes. 
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Detectability at 1.9 lambda/D vs used images number 
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Figure 3. SCC (plus and diamonds) and differential imaging (crosses and trian- 
gles) detectabilities at 5a at 1-9-^ versus the number of used images for Inm 
andbnm rms non-common static aberrations (respec. dotted and dashed line). Pho- 
ton noise (solid line) at 1.9-^ is also plotted. Log-log scale is used. The SCC reference 
image exposure time is 18 s. 

5.2 Detectability versus number of used images 



To check that the SCC noise decreases as the square root of the number of 
images recorded (following the photon noise), we plotted on figure 3 the SCC 
detectability (plus and diamonds) at 1.9-^ versus the number of used images. 
In comparison, we added the same plot for the differential imaging detectabil- 
ity (crosses and triangles). We also plotted theoretical photon noise (solid line) 
and non-common static aberrations noises (dotted and dashed lines for respec- 
tively 1 nm and 5nm rms). As seen in figure 2, the differential imaging, unlike 
the SCC, is limited by non common static aberrations. Furthermore, where Ir 
is greater than one photon, the SCC image is limited by photon noise and 
decreases as the square root of the number of used images. However, there is a 
minimum detectability that the Self-Coherent Camera can reach. Indeed, even 
if the photon noise is minimized on 1^ using an 18 s exposure time, it is still 
present and limits the SCC detectability (flat evolution after 5000 images). 
The exact impact of the reference beam noise is still under study. 
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Seff^Cah^mnt Cftmersr IJlllkt^mAl imaging 




Figure 4. Final Self-Coherent Camera (left) and differential (right) images using a 
single (top) and 400 images (bottom). The companion, l-9-^ separated, is 10^ less 
bright than the host star. The SCC reference image exposure time is 18s. An algo- 
rithm based on cosmic removal is used to correct the SCC images from spikes due 
to the division by Ir. Power scale is the same for the four images. 



5.3 Resulting images 



Prom figure 4, it looks like 400 SCC images (corresponding to only 2.4 s) 

are sufficient to detect a 10~^ companion at an angular separation of 1.9-^ 
from the host star. This is true whatever the static aberration amplitude is. 
On the contrary, the differential imaging technique is limited by non-common 
static aberrations. We simulated a companion-star system and observed the 
corresponding images (figure 4). As expected, the companion is not detected 
on the final images using a single image for both techniques that are photon 
noise limited (top of the figure). Using 400 images, the companion is lost in 
non-common static speckles on the final differential image whereas it is well 
detected on the final SCC image. 



9 



6 Conclusion 



In this paper, we tiave recalled the principle of the Self-Coherent Camera 
(SCC). We have described the numerical simulation we developed to eval- 
uate the performance of the SCC. We found that, as opposed to standard 
differential imaging, the SCC is not limited by the aberrations in the optical 
elements. The detection limit with the SCC roughly follow the photon noise 
limitation (cf. figure 2 and 3). Thus, to detect a faint companion very close 
to its host star (10~^ at 1.9-^), only a short observation time is mandatory 
(about 2.4 s for a magnitude 5 star on an 8- meter telescope). To detect fainter 
companion as a 10~^° exoplanet, the coupling of coronagraphs and the SCC 
will be necessary. This coupling is under study. We are also analyzing differ- 
ent data processing technique to minimize the impact of the division by the 
reference image Ir. A laboratory experiment is also under development to 
compare the expected capabilities presented in this paper with the effective 
performance. 

References 

[I] Baudoz, P., Boccaletti, A., Baudrand, J., Rouan, D., 2006, Proc. of lAU 
Colloquium 200, 553 

[2] Schneider J., The Extrasolar Planets Encyclopaedia, http://exoplanet.eu/ 

[3] Baraffe, I., Chabrier, G., Barman, T. S., Allard, F., Hauschildt, P. H. 2003, 
A&A, 402, 701 

[4] Guyon, O., Pluzhnik, E., Kuchner, M., Collins, B., Ridgway, S., 2006, ApJ Sup., 
167, 81 

[5] Marois, C., Doyon, R., Nadeau, D., Racine, R., Riopel, M. Vallee, P., Lafreniere, 
D. 2005, PASP, 117, 745 

[6] Racine, R., Walker, G., Nadeau, D., Doyon, R., Marois, C. 1999, PASP, 111, 
587 

[7] Seager, S., Whitney, B., Sasselov, S. 2000, ApJ, 540, 504 
[8] Baba, N., Murakami, N. 2003, PASP, 115, 1363 
[9] Codona, J., Angel, R. 2004, ApJ, 604, L117 
[10] Guyon, O. 2004, ApJ, 615, 562 

[II] Labeyrie, A. 2004, EAS Conf. Puh. Series, 12, 3 

[12] Ollivier, M., Mariotti, J.-M. 1997, Appl. Opt., 36, 5340 



10 



[13] Wynne C. 1979, Opt. Comm., 28, 21 

[14] M. Sarazin, A. Tokovinin, 2002, ESO Conf. and Workshop Proc, 58, 321 

[15] F. Rigaut, J.-P. Veran, O. Lai, 1998, SPIE Proc, 3353, 1038 

[16] L. Jolissaint, J.P. Veran, 2002, ESO Conf. and Workshop Proc, 58, 201 

[17] C. Cavarroc, A. Boccaletti, P. Baudoz, T. Fusco, and D. Rouan, 2006, A&A, 447, 
397 

[18] A. Duparre, J. Ferre-Borrull, S. Gliech et al., 2002 Appl. Opt., 41, 154 



11 



